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Summary 

This paper presents summary details of an aperiodic time asymptotic numer- 
ical solution for the square drive cavity at Re = 10000. The data presented is for 
6100 < t < 7100, and is representative of the data that characterizes the aperi- 
odic asymptotic state. Complete details are not given in this paper for either the 
transient evolution or the asymptotic state. 
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Introduction 

A persisting final oscillation in the aspect ratio two driven cavity was shown 
in [5] at Re — 10000 and t ss 300. A periodic time asymptotic flow in the aspect 
ratio two driven cavity was presented in [3] at Re = 5000 and t « 3800. This 
paper presents summary details of an aperiodic time asymptotic numerical solution 
for the square drive cavity at Re = 10000 and t « 6000. Further details are in 
[2]. Unsteady time asymptotic calculations for the square regularized driven cavity 
have been presented in [7]. 


Numerical Method 

The velocity field for two dimensional time dependent incompressible fluid flow 
in a spatial domain H may be written in terms of the streamfunction tp{x,t) as 
u(x,t) = and v (x,f) = The flow equations may be written as 


d Atp 

dt 


1 2 dtp dtp dtp 

tf ’ + "di^lh, ~ ~fy A 


dtp_ 

dx' 


for x in O, and t > 0. 


( 1 ) 


The data for the impulsively started driven cavity consists of the initial values 
0(x, 0) = 0 for x in fl and t = 0, and the boundary conditions tp(x,t) = 0 for x in 
dQ and t > 0, and f^(x, t) = 1 for x on the cavity lid and t > 0, and |*-(x, t) = 0 

for x on the cavity walls and t > 0, where is differentiation in the exterior normal 
direction at the boundary, and where the lid is moving from left to right for t > 0. 

Let a discrete approximation for the streamfunction on a square uniform grid 
be z = {z m : m = 0, 1, . . . }. If La is the conventional five point centered difference 
approximation to the Laplacian, if Bi is the conventional thirteen point centered 
difference approximation to the Biharmonic operator, and if 6 X and 8 y are the 
conventional centered difference operators, then a discretization of equation (1) is 


_ At , 3A t 

_La(i 2i2« Bl Z 2~ 


i, ^,(r)La(i")j - 6, ^.(i")La(£')j 


( 2 ) 




This algorithm is second order accurate in both time and space. Standard finite 
difference approximations are used to incorporate the boundary condition. The 
velocity components are directly recovered from the discrete streamfunction solution 
using central difference approximations for the and The velocity solution 
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is exactly discretely divergence free with respect to a central difference formulation 
of the mass conservation equation. Further details are given in [2] and [4]. 

The solution at time t n+x of the implicit equation (2) is obtained by a multigrid 
method, as in [2], The Biharmonic operator is factored as two Laplacians, point 
Gauss-Seidel relaxation is used for the smoothing operator, and linear restriction 
and prolongation operators are used. A V-cycle iteration scheme is used, with 3 
iterations per grid level while coarsening, and none while refining. At each time 
step, 10 to 15 iteration cycles are used to reduce residuals to less than 5.0 x 10" 12 . 

Numerical Results 

Computations are reported for the impulsively started square driven cavity at 
Re = 10000 on a uniform 128 x 128 grid with At = These computations 

show that the time asymptotic state for this combination of algorithm, data, and 
grid is an aperiodic laminar flow. Among the indicators that have been used [3] 
for tracking convergence to an asymptotic state are the relative L x norm of the 
streamfunction change per time step / 1 — V'"/ \^ij the total ki- 

netic energy |AxAy£\ . ||u t n + 1 1| 2 , and the minimum streamfunction value, each 
of which is obtained from the data on the interior grid. 

The relative L x norm for the streamfunction change per time step is 0(10 6 ) 
for t > 100, so that virtually the entire time evolution of the system and its time 
asymptotic state are in a range of very small change per time step with respect to 
the streamfunction solution surface as a whole. The initial transients completely 
dominate the time evolution of this system until at least t ta 500. The relative L x 
norm for the streamfunction change per time step seems to approach its asymptotic 
pattern by t « 1600, with persisting variation in the range [0.00006,0.00009]. The 
total kinetic energy of the flow on the interior grid in the cavity is greater than 
0.036 by t « 600, but there is clearly evident transient change until at least t « 1900 
when the asymptotic evolutionry pattern is approached. Persisting variation for the 
total kinetic energy is in the range [0.03624,0.03630], The minimum streamfunction 
value approaches its asymptotic time evolution pattern and range by t « 1900, with 
persisting variation in the range [—0.10823,-0.10804]. The global streamfunction 
minimum is also point data, since the global minimum for the asymptotic state is 
at a fixed point in the very center of the primary circulation. The combination 
these three indicators suggests that the asymptotic attractor for this flow has been 
approached by 1 « 2000. 

There were small transient changes in phase portraits for this data, so the 
computation was run until t = 7100 in order to ensure that the asymptotic state 
had indeed been reached. There is a period of T « 200 for the asymptotic flow. 
The resolution of the asymptotic attractor with phase portraits and spectral analysis 
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required at least 1000 nondimensional time units, or data from 160000 time steps. 
Figure 1 shows the time series data for the relative L x norm of the streamfunction 
change per time step, and Figure 2 shows the streamfunction minimum, both for 
6100 < t < 7100. This data shows approximately five complete periods of the 
asymptotic state. Notice the complicated beating of multiple frequencies. Figure 3 
shows the power spectrum for the total kinetic energy, for normalized frequencies 
less than 5. This spectrum is representative of the complicated spectra for the data 
from this flow, with multiple discrete spectral lines, with a primary frequency of 
c p « 0.237, with various harmonics, and with a minimal significant frequency of 
c m ~ 0.005. For the relative L t norm of the streamfunction change per time step 
the ratio of spectral weights for c p and c m is 61.3, for the streamfunction minimum 
it is 132.0, and for the total kinetic energy it is 5.66. It seems likely that there 
are at least two incomensurate frequencies with c m approximately equal to their 
difference. The power spectra for all three of the scalar time series also reflect 
a complicated substructure in their spectral lines, with the peaks in the power 
spectrum away from c p tending to spread into a cluster of subpeaks, with frequency 
seperations that are approxunately 0.005. Figure 4 shows a phase portrait of the 
relative Li norm of the streamfunction change per time step versus the total kinetic 
energy, for 6100 < t < 7100. This time interval was required to define the entire 
phase portrait. This phase portrait clearly shows a complex but regular dynamical 
process. Figure 5 shows three representative streamfunction contour plots for the 
asymptotic flow . The unsteady dynamics are laminar, and are visibly concentrated 
in the relatively weak secondary flow structures. In each of the secondary flow 
structures weak tertiary circulations form underneath the separated primary flow 
near the point of separation, and then are convected into the secondary flow as a 
whole, to interact with and join the focus of the secondary flow. 

Discussion 

There are well known steady solutions for the square driven cavity at Re = 
10000 in [1] and [6]. These steady results do not contradict the current work since 
they were obtained with steady state codes. Theoretical expectations [8], the un- 
steady asymptotic results in [7] for the square regularized driven cavity using a 
Chebyshev-Tau approximation, and the results in [3] for the rectangular driven cav- 
ity using the same method as in this work all support these unsteady asymptotic 
results in the square driven cavity at Re = 10000. Continuing work in preparation 
[2] suggests that the flow in a square driven cavity has a Hopf bifurcation to a 
periodic asymptotic state for Re < 9000. 
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